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The recent observations of the massive pulsars PSR J1614-2230 and of PSR J0348+0432 with 
about two solar masses imply strong constraints on the properties of dense matter in the core of 
compact stars. Effective models of QCD aiming to describe neutron star matter can thereby be 
considerably constrained. In this context, a chiral quark-meson model based on a SU(3) linear 
(T-model with a vacuum pressure and vector meson exchange is discussed in this work. The impact 
of its various terms and parameters on the equation of state and the maximum mass of compact 
stars are delineated to check whether pure quark stars with two solar masses are feasible within this 
approach. Large vector meson coupling constant and a small vacuum pressure allow for maximum 
masses of two or more solar masses. However, pure quark stars made of absolutely stable strange 
quark matter, so called strange stars, turn out to be restricted to a quite small parameter range. 


I. INTRODUCTION 

The core of a massive star which has exhausted its 
thermonuclear fuel collapses in a supernova which is one 
of the most extreme events known to occur in the uni¬ 
verse. The relic in the aftermath of such a cataclysmic 
event constitutes also one of the most extreme objects 
known to exist in the universe. Being more massive than 
our sun, but only around 20 to 30 kilometers in diameter, 
these so called neutron stars harbor the densest material 
known at present. The density in the core of such an 
object can even exceed nuclear density (po ~ 2.5 • 10^^ 
g/cm^). 

For many decades it has been speculated that the ma¬ 
terial in the core of neutron stars might consist of quark 
matter 0,0. Neutron stars, or better compact stars, 
with a core consisting of quark matter are dubbed hy¬ 
brid stars. Compact stars which are entirely made of 
quark matter, besides maybe a small layer of a crust of 
nuclei, are so called strange stars 00. Strange stars 
can only be realized in nature if strange quark matter 
is absolutely stable, i.e. the true ground state of matter 
00. Hybrid stars as well as strange stars have been 
usually modeled by using the MIT bag model, see 0 
for a review, also with corrections from hard-dense-loop 
calculations 0. The Nambu-Jona-Lasinio (NJL) model 
has been adopted for describing quark matter in compact 
stars by using scalar meson helds only in 0110. Vector 
meson helds have been added within the NJL-model in 
[HI showing that their contribution can substantially in¬ 
crease the maximum mass of a compact star with quark 
matter. Effects from color-superconductivity have been 
also considered for the properties of compact stars within 
NJL-type models, see [Q for a review. We point out that 
it was well known that hybrid stars can be as massive as 
two solar masses in various approaches [3 [a and could 
even masquerade as neutron stars [Tsj l. 

The recent measurements of the masses of the pulsars 
PSR J1614-2230 [10 and of PSR J0348-k0432 [l|] with 


M = 2.01 ± 0.04Mq impose considerable constraints on 
the equation of state (EoS) for compact stars. Pure quark 
stars based on the simple MIT bag model could be ruled 
out on the basis of this measurement unless additional 
terms from an effective one-gluon exchange or from color- 
superconductivity are incorporated Hip. However, 
the MIT bag model fails in describing QCD lattice data 
at nonvanishing temperature questioning its applicability 
for describing dense quark matter [20. Effective models 
of QCD based on chiral symmetry are able to describe 
the lattice data at nonvanishing temperatures, as the 
Polyakov-loop extended versions of the NJL-model 1^ 
or the Polyakov-loop quark-meson (PQM) model (22 h^ . 
NJL-type models for hybrid stars have been investigated 
for being compatible with a maximum mass of com pac t 
stars of at least two solar masses by several groups (2^ 
[10. Recently, the connection between NJL-type models, 
the Schwinger-Dyson approach and an extended version 
of the MIT bag model with vector-like interaction terms 
has been pointed out in [33 - [30 for modeling compact 
star matter incorporating features of QCD. To our knowl¬ 
edge, the quark-meson model has not been used to inves¬ 
tigate the properties of hybrid stars or quark stars so far, 
in particular not in view of the now well established two 
solar mass limit for compact stars. 

In this work, a modified linear-cr-model is used to de¬ 
scribe compact star matter consisting of quark matter 
only. The linear cr-model is well suited to consider the 
chiral symmetry breaking patterns of strong interactions 
[30 . The quark-quark interaction is mediated by the ex¬ 
change of meson fields. These interactions are conceptu¬ 
ally different from the NJL model which considers point 
coupling terms between quarks. We consider scalar- and 
vector meson contributions to effectively model the at¬ 
tractive and repulsive character of the strong interac¬ 
tion. The maximum masses of pure quark stars are 
calculated by solving the Tolman-Oppenheimer-Volkoff 
(TOV) equations for different choices of the parameters 
of the quark-meson model. We further investigate the 
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parameter space required for the existence of absolutely 
stable strange quark matter and strange stars and con¬ 
front these results with the new maximum mass limit 
for compact stars. This model has been studied in [s^ 
for the properties of quark matter under the conditions 
present in core-collapse supernovae. However, the two 
solar mass constraint has not been investigated so far. 

II. THE CHIRAL QUARK-MESON MODEL 


The indices n and s denote the nonstrange and strange 
quark contributions. The quark fields couple to the 
scalar- and vector meson fields (t„, as, uj, p, and (j) with 
the respective coupling strength pi, which are related by 
SU(3) flavor symmetry to one overall coupling constant 
for the scalar meson and to another one for the vector 
coupling constant gy. 

Considering stationary fields in mean field approxi¬ 
mation, the derivative terms in Cm vanish and the La- 
grangian is given by 


The quark-meson model is based on the linear cr-model 
as discussed in detail in [13, and couples mesons and 
quarks by utilizing chiral symmetry. The mesonic con¬ 
tribution in this model is given by 

Cm = tr{d^ip)'^{d^ip) + 

- \i[tr{Lp'^p)Y - \2tr{ip^ipY 

— (triV^V)) 

- tr[H{'p + (^'1')] -I- c (det(i^'*') -I- det((/5)) (1) 

for SU{3) X SU{3) chiral symmetry incorporating the 
scalar {(p) and vector (V^) meson nonet. Here, stands 
for the vacuum mass of the vector mesons and Ai, A 2 , tuq, 
and c are the standard parameters of the linear cr model 
to be fixed below. The matrix H describes the explicit 
breaking of chiral symmetry. The quarks couple to the 
meson fields via Yukawa-type interaction terms 

Cq = ^ {i0 - g,pip - gvl^Vf,) (2) 

with the couplings strengths g^ and pv for scalar and vec¬ 
tor mesons, respectively. Both contributions are forming 
the SU(3) Lagrangian C = Cm + Cq of the chiral quark- 
meson model. 

In the mean-field approximation, the matrix tp consists 
of just the scalar nonstrange field and the strange field 
<7s'- 


C = Cm + Cq 

= \ + mlp^ + rnlcj)'^) 

- ^(O'n + <^^8? - ^(o-n + 

2 

- ^(o'n + 0-s) + V2alasC + h„an + hsUs - B 

+ {i0 - guj"f°UJ - PpTJ^p - Pnan) 5'n 

+ ^s{i0-PsiJs- ( 6 ) 


Here, a vacuum energy term B has been introduced in 
addition, see the discussion in [3, |4l|, . The electrons 

will be treated as a free noninteracting Fermi-gas. The 
potential of the Lagrangian ([5)) then reads 


V = - i + mlp^ + 

+ ^(^n + + '^(^n + 0's) 

+ + ^s) - • C - hnCTn “ ^s^s + B (7) 


for the meson fields. 


HI. THE EQUATION OF STATE 



p = 

and the vector fields are described by 

V2 



(3) 


(4) 


where p stands for the zeroth component of the isovector 
vector field, uj for the nonstrange vector field and tp for 
the strange vector field, assuming ideal mixing. Note 
that the spatial components of vanish in the mean- 
field approximation in the static case to be considered 
here[^. The Lagrangian for the quark fields then reads: 

{i0 - - PpTJ^p - 5„CT„) 

+{i0 - PsCTs - g.pj^pj'^s (5) 


The grand canonical potential is related to the parti¬ 
tion function via 


n = - 


InZ 

~T 


= -p- 


( 8 ) 


The partition function Z can be computed by a Feynman 
path integral over the quark fields. Performing the inte¬ 
gration in mean field approximation, the classical station¬ 
ary mesonic background Helds can be replaced by their 
non-vanishing vacuum expectation values: 


Z = / n BaaBwa I U ^"^£+'*'7%'!')) 

a 

and one arrives at 


Q = V - ^ / k^dk ■ M 

Jo 


(9) 

( 10 ) 




3 


where V is the potential given in eq. 0 and the short¬ 
hand notation for A4 is 


/ —Bu+y-u—gu!'^ — gpp 

M=\n{l + e - "F - 


-I- In ( 1 -I- e 


-Ed + Pd~9^^ + 9pP 


).!„( 

)+!„( 


— Eu — pu-\-guj^-hgpp '' 

iH-e-- 


1 + e' 


~Ed-Pd+9i^^-gpp ' 


~^s'^Ps~gcj}<f> 

-f In 1 -h e T 


— Es—ps+g^<p\ 

+ In 1 -f e-3^- 


( 11 ) 


where the flavor dependent energy is 


the isovector number densities and the strange number 
densities due to SU(3) flavor symmetry. 

The coupled equations of motion of the meson fields 
have to be solved self-consistently. Using eq. ([5]) and the 
relation 

e = U - ^ PiHi (22) 

i 

the resulting field values then determine the EoS, which 
serves as an input to solve the TOV-equations: 


Ef = \Jk} + {9f(FfY (12) 

Compact star matter can be treated in the zero temper¬ 
ature limit T —>■ 0. The equations of motion of the five 
meson fields are given by minimizing the thermodynamic 
potential: 


dm 

dr 

dp 

dr 


47rr^e(r) 


G€(r)m(r) 

{cry 


X 




47rr^p(r) 

m{r)c^ 


/ 2Gm{r) 

\ c^r 


(23) 

(24) 

(25) 


dn _ an _ dn _ dn _ on 

dan das doj dp dcj) 

For the scalar mesons one finds the gap-equations 


(13) 


for the mass radius relation of compact stars. 


A. Parameter range 


^ = \ian{al + a^) + ycr^ -f TO^„cr„ - hn (14) 

3^ f dk-e dk-e \ ^ 

\J^ En Jo Ea ) 


Within the mean field approximation the mesonic fields 
will be treated as classical background fields. In the vac¬ 
uum they will be replaced by their vacuum expectation 
values (VEV). In the scalar sector they are determined 
by the weak decay constants 


and 


dVt 

das 


= Xias{ai + ai) + X2ai+mlas-hs (15) 


iglas f dk ■ k^ 


Es 


= 0 


and for the vector fields 


dn 

duj 

dfl 


= -rris^uJ + 




2 I ^9p 


I dk ■ J dk ■ k^ j = 0 (16) 

pkp pkp \ 

' dk-k"^- dk-k^] = 0 (17) 


dn 

dcj) 


= 




dfc • = 0 


(18) 


and the respective Fermi momenta are 


ku = ^ - gpP^ - (ffnCr„)2 (jg) 

kd = y (m, + y - + gpP^ - {9n(FnY (20) 

ks = y (m, + y - g^ 4 >) - {gsCTsY (21) 


to guarantee charge neutrality through the electron 
chemical potential Pe- Note that the source terms for 
the vector fields are given by the vector number densities, 


< an >= U = 92.4 MeV (26) 

< as >= ~ = 94.47 MeV (27) 

-v/2 

with Jk = 159.8/-\/2 MeV. The model incorporates four 
free parameters. The constituent quark mass will be var¬ 
ied in a range around = 100 MeV to 400 MeV. It de¬ 
termines the scalar coupling for the nonstrange (E51) and 
strange (HID condensates via the Goldberger-Treiman- 
relation 

m„ 

gu,d = 9n = (28) 

/tt 

and from SU(3) symmetry 

5s = • a/2 . (29) 

The vector coupling is independent on the constituent 
quark mass, it will be varied in a scale similar to the one 
of the scalar coupling, ~ pn, to study its influences in 
an appropriate range |38l| . The coupling constant of the 
((>-meson is fixed again by SU(3) symmetry: 


5.=5p = ^ (30) 

To fix the remaining parameters of the scalar mesons, the 
following masses have been adopted: 

= 138 MeV (31) 

mK = 496 MeV (32) 

mn = 547.5 MeV (33) 

= 957.78 MeV (34) 
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The experimentally not well determined mass of the cr- 
meson could cover a range from rria = 400 MeV to 
1000 MeV. The explicit symmetry breaking terms are 
defined via the Gell-Mann-Oakes-Renner relation [s^ as 
follows 

hn = Uml (35) 

hs = (36) 

The parameter Ai has to be determined via the relation 


Ai) O- m^(m^(Ai), Ai) o m^(Ai) (37) 

and A 2 is given by 

^ ^ 3m|,(2/j^: - /^) - m^(2/x + U) 

" m + SfKifK-fn))ifK-f.) 

m + ifK{fK-U))UK-h) ^ ’ 


The large mass of the r/-meson (which, as a Goldstone 
boson should be nearly massless, see (^,13) t)e im¬ 
plemented by the axial anomaly term with the coupling 
strength c determined by 


c = 


2 2 

Ik - U 


A 2 • (2/^ - Jtv)- 


(39) 


The last free parameter is the vacuum pressure at van¬ 
ishing chemical potential B, which we choose to cover 
a range from = 0 MeV to 140 MeV. While vary¬ 
ing one parameter in the following, the other parame¬ 
ters will be held fixed at some canonical values chosen 
to be rricr = 600 MeV, ruq = 300 MeV, = 2 and 
51/4 = 120 MeV. 



FIG. 1. The nonstrange scalar condensate for various values 
of the vector coupling constant g^i at fixed values of iriq = 
300 MeV, = 600 MeV and = 120 MeV. 



IV. RESULTS 


FIG. 2. The strange scalar condensates for various values 
of the vector coupling constant at fixed niq = 300 MeV, 
nicr = 600 MeV and = 120 MeV. 


A. Variation of the vector coupling constant 

The vector coupling models the repulsive character of 
the strong interaction. Analogue to the scalar coupling, 
which is dependent on the constituent quark mass, the 
vector coupling will be changed in a range from = 1 
to 7. 

Figure [T] shows the nonstrange scalar condensates for 
various values of g^j. There is a crossover for values of 
51 ^ ^ 4. For values of ^ 4 there is a jump in the 
scalar condensate, which is present at smaller chemical 
potential fiq for smaller values of guj. The chiral transi¬ 
tion for g^; = 2.0 takes place at /r ^ 325 MeV. Because 
of explicit symmetry breaking, the condensate does not 
vanish entirel y in the chirally restored phase, see e.g. the 
discussion in [37|. 

Note that the jump is also visible for the strange con¬ 
densate, which can be seen in Fig. [51 Because of the 
larger mass of the strange quark the jump is not that 


pronounced and the field changes rather smoothly stay¬ 
ing at larger values even in the chirally restored phase up 
to the largest considered values of the chemical potential. 

An increase of the repulsive coupling g^j provokes the 
scalar fields to increase too. This is due to the sub¬ 
struction of the terms including the vector fields from 
the chemical potential /ig in eqs. (1191) . (1501) and (1211) . 
Larger vector field terms require larger chemical poten¬ 
tial to compensate for their impact. 

The EoS for various values g^j is plotted in Fig.|3l For a 
larger value of the vector coupling constant, the repulsion 
between quarks increases, the slope of the EoS increases 
and the EoS becomes stiffer, i.e. quark matter requires 
more pressure to be compressed to a given energy density. 

The corresponding solutions of the TOV equations can 
be seen in Fig. 01 With a stiffer EoS the maximum mass 
increases. This increase in the maximum mass is due to 
the repulsive character of the vector mesons which gives 
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FIG. 3. The EoS for various values of the vector coupling 
constant g^i at fixed mq = 300 MeV, nia = 600 MeV and 
5^/^ = 120 MeV. 



FIG. 4. The mass-radius relation for various values of the 
vector coupling constant g^j at fixed rriq = 300 MeV, = 
600 MeV and = 120 MeV. 

a higher pressure at a given energy density (see Fig. [3]) 
and therefore quark matter is able to stabilize more mass 
against the pull of gravity. With the choice of gqj = 4.0 
and the other parameters held fixed at ruq = 300 MeV, 
rua = 600 MeV and = 120 MeV one reaches a max¬ 
imum mass of ^ 2.1Mq at a radius of 10 km. Smaller 
values of the vector coupling constant result in compact 
star configurations with a maximum mass smaller than 
2Mq which is in conflict with observations. Hence, the 
vector-like interactions between quarks are necessary in 
the model used to achieve the 2Mq mass limit in a phys¬ 
ically reasonable range of our parameters. 

B. Variation of the constituent quark mass 

Raising the value of the constituent quark-mass in¬ 
creases the scalar coupling, see equations (1351) and (1351) . 



FIG. 5. The EoS while varying mq with fixed = 2.0, 
mcr = 600 MeV and = 120 MeV. 



R[km] 


FIG. 6. The mass-radius relation for various values of mq at 
fixed guj = 2.0, ma = 600 MeV and B^^'^ — 120 MeV 


The values of the scalar condensate have been checked 
with the results of and found to be in accordance. 

We find that a crossover transition is present for values 
of fJ,q ^ 300 MeV, whereas for values fiq > 300 MeV a first 
order phase transition emerges, and the chiral condensate 
jumps. For even higher values, the strength of the first 
order chiral phase transition increases. 

The EoS shown in Fig. [5] displays a softening for 
increasing values of mq, leading to smaller maximum 
masses of the compact star configuration, as can be 
seen in Fig. [51 The trend of the curve of the EoS for 
mq > 300 MeV shows a slightly different behaviour. For 
mq = 100 MeV the 2Mq -limit is exceeded. In this case 
the scalar condensates would exhibit a smooth crossover¬ 
like behavior as a function of the chemical potential. 
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FIG. 7. The EoS for various values of mo- at fixed gu; = 2.0, 
niq — 300 MeV and = 120 MeV. 



R[km] 


FIG. 8. The mass-radius relation for various values of mo at 
fixed gu, = 2.0, m^ = 300 MeV and = 120 MeV. 


C. Variation of the u-Meson mass 

We find that with increasing a meson mass rria the 
phase transition becomes a crossover (mo > 800 MeV), 
whereas values of < 800 MeV lead to first order phase 
transitions at /r, ~ 300 MeV. For mo = 600 MeV the first 
order phase transition takes place at ^ 330 MeV. The 
behavior of the nonstrange and strange condensates is 
found to be similar to the discussion in section IIV Al 

The resulting EoS is shown in Fig. [T] The EoS softens 
with increasing scalar meson mass m„. Hence, one ex¬ 
pects a mass-radius relation which is located at smaller 
values of the mass and radius for increasing the a me¬ 
son mass, which can clearly be observed in Fig. [S) So a 
smaller value of leads to higher maximum masses of 
the compact star. 

The mass of the sigma meson m„ is directly related to 
the parameters Ai and A 2 , according to (1571) and (1551) . 
These parameters are, among others in the SU(3) case, 



FIG. 9. The equation of state for various values of the vacuum 
constant B^^'^ at fixed values of nia — 600 MeV, g^, = 2.0 and 
niq — 300 MeV. 


mainly responsible for the potential depth for sponta¬ 
neous symmetry breaking. The impact of the sigma me¬ 
son mass on the chiral condensate and the EoS is nontriv¬ 
ial. A simple explanation might be that for a higher mass 
of the scalar meson more energy is needed to overcome 
the barrier and reach the second minimum of the poten¬ 
tial, which leads also to a shift of the chiral phase transi¬ 
tion to a larger chemical potential fig and a crossover-like 
behavior of the chiral condensates at high densities. 


D. Variation of the vacuum constant 

Figure O shows the equation of state for various val¬ 
ues of the vacuum constant The equation of state 

does not change significantly when varying the vacuum 
constant. The vacuum constant B drops out in the equa¬ 
tions of motion (1751) so that it does not affect the values 
of the meson fields. 

The slopes of the curves are very similar to each other. 
The pressure vanishes at lower energy densities for a 
smaller vacuum constant. Note that the curve for the 
vacuum constant B^^'^ = 40 MeV is not simply shifted 
but the jump in the energy density is in this case at a 
small, but non-vanishing value of the pressure. This can 
be observed in the inlay of Fig. [51 This property leads to 
an additional stable branch in the mass-radius relation 
shown in Fig. [TO] For = 40 MeV there are two sta¬ 
ble branches. One up to 0.6 Mq at a radius of ~ 19 km 
and a second branch between ~ 14 km and ~ 12 km ra¬ 
dius and up to a maximum mass of ^ 2Mq. These so 
called twin star solutions Bin are beyond the scope of 
this article and are discussed in more detail in [8|. 
and in the forthcoming publication [s^. See (45j| for a 
detained discussion on stability. 

The mass-radius relation depicted in Fig. [TU] demon¬ 
strates that larger maximum masses can be achieved for 
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R[km] 



FIG. 10. The mass-radius relation for various values of the 
vacuum constant with fixed values of m,, = OOOMeV, 

gu) = 2.0 and rriq = 300 MeV. The maximum mass increases 
significantly for a lower value of the vacuum constant B. 

smaller values of the vacuum constant B. We point 
out that the standard choice of the vacuum pressure in 
refs. [1, Hlj was a value of = 120 MeV. For this 
value the resulting mass-radius relation has a maximum 
mass around I. 6 M 0 at a radius of around 9 km, when 
fixing the other parameters of the model at their stan¬ 
dard values. The vacuum constant has a strong influence 
on the mass-radius relation, and the maximum mass of a 
compact star can easily reach values of 2Mq. Choosing 
B^/^ = 40 MeV the resulting compact star configurations 
have a maximum mass of ^ 2Mq at a radius of 12 km. 

E. Stability considerations 

The properties of quark matter have to fulfill certain 
conditions in order to allow for the possibility of pure 
quark matter stars studied above. 

Normal matter, consisting of ordinary nuclei, is sta¬ 
ble on cosmological timescales, so it does not decay to 
quark matter with its quark constituents of up-quarks 
and down-quarks. This observation requires that two fla¬ 
vor quark matter can not be more stable than ordinary 
nuclear matter, meaning that the energy per baryon has 
to be higher than the one of the most stable known el¬ 
ement in nature: ®®Fe. We adopt a value of energy per 
baryon of E/A = 930 MeV for nuclei and add a 4 MeV 
correction due to surface effects of lumps of quark matter 
balls as discussed in . Hence the critical condition for 
two-flavor quark matter reads 

>311 MeV (40) 

p=0 

This condition guarantees the stability of atomic nuclei, 
meaning that atomic nuclei do not dissolve into their con¬ 
stituent quarks. 


FIG. 11. Maximum masses for three-flavor pure quark stars in 
the plane of the vacuum constant B and the vector coupling 
constant The two-flavor and three-flavor lines delineate 
the region for hybrid stars (above the two-flavor line) and 
pure quark stars (below the three-flavor line). Pure quark 
star configurations do not appear for the parameters chosen 
(rrio- = 600 MeV and niq = 300 MeV). 

On the other hand three-flavor quark matter, i.e. quark 
matter consisting of up-, down- and strange quarks, could 
be more stable than ordinary nuclei which is the Bodmer- 
Witten hypothesis Ordinary nuclear matter can 

not decay to strange quark matter, as there is a barrier 
between these two states of matter due to strangeness 
conversion via weak interactions with a corresponding 
conversion timescale much longer than the age of the 
universe. The presence of the new degree of freedom, 
the strange quark, in quark matter lowers the overall en¬ 
ergy per baryon, so that this state could be energeti¬ 
cally favorable compared to nuclear matter. Hence, the 
condition for absolutely stable strange quark matter of 
E/A < 930 MeV can be recast in the form 

< 310 MeV (41) 

p=0 

In the following we will denote the physical condition for 
two-flavor quark matter, eq. dMl), the two-flavor condi¬ 
tion (or two-flavor line in the plots) and the one for three- 
flavor quark matter, eq. m, the three-flavor condition 
(or three-flavor line in the plots). 

Figures ITlI and fT^ depict the maximum masses in de¬ 
pendence of the parameters 5^/^, mo- and g^j. In general, 
smaller values of the vacuum constant B and higher val¬ 
ues for the vector coupling constant g^j lead to higher 
maximum masses. 

Figure [TT] shows as a contour plot the dependencies 
of the maximum mass of pure quark star configurations 
on the vacuum constant B'^/^ and the vector coupling 
constant g^j. For B^/'^ = 120 MeV and g^, = 2.0 one 
finds a maximum mass of about I. 6 M 0 , which can be 
cross-checked with Figs, m IHl [5] and [TUI 
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FIG. 12. Maximum masses for three-flavor pure quark stars in 
the plane of the vacuum constant and the sigma meson 
mass rria. The two-flavor and three-flavor lines delineate the 
region for hybrid stars (above the two-flavor line) and pure 
quark stars (below the three-flavor line). Pure quark star 
configurations are possible for a small window at large values 
of the sigma meson mass rria (here = 2 and niq = 300 MeV 
are held fixed). 


The two-flavor line indicates the two-flavor constraint 
(TO and the three-flavor line the three-flavor constraint 
(IdlT) . In the area above the two-flavor line the condi¬ 
tion dSl) is fulfilled, i.e. normal matter can not decay 
to two-flavor quark matter as observed in nature. For 
a vanishing vector repulsion the vacuum constant has to 
be larger than > 60 MeV in order for quark matter 
to obey the two-flavor constraint. The critical value for 
decreases towards = 20 MeV with increasing 
vector coupling constant. Hybrid star configurations are 
located in the parameter range above the two-flavor line. 

In the area below the three-flavor line the condition 
(l4T|l is fulfilled. Within the chosen parameters the three- 
flavor line is nearly independent on the vector coupling 
constant g^j being between B^^'^ ^ 28 MeV for a van¬ 
ishing vector coupling constant and B^^^ ~ 20 MeV for 
g^j = 9.0. The vacuum pressure dictates at which en¬ 
ergy density eo the pressure vanishes, which determines 
via the Hugenholtz-van Hove theorem the binding en¬ 
ergy of quark matter: Eb/A = = ^olnB- As, there is 

no region, where both conditions (HOI) and dm are ful¬ 
filled simultaneously, no pure quark star configurations 
emerge in the contour plot for the chosen parameters of 
TOct = 600 MeV and = 300 MeV. 

Figure \V3 shows the maximum masses and the quark 
matter constraints (l40l) and (IdTI) in the parameter plane 
of the vacuum constant and the cr-meson mass 

Too-. The two-flavor line starts at B^^^ = 115 MeV for 
rria = 400 MeV, with the three-flavor line being slightly 
below the two-flavor line. Both lines decrease in a sim¬ 
ilar manner with increasing mass of the cr-meson. At 
B^/^ = 26 MeV and = 725 MeV both lines inter¬ 
sect and from that point on the three-flavor line is above 


FIG. 13. Maximum masses and stability configurations for 
pure quark star configurations as in Fig.[T2]for a slightly larger 
vector coupling constant of g^i = 3.0. 


the two-flavor line. This means that from there on, in a 
small parameter space (the small area enclosed by both 
lines), pure quark star configurations are stable. The cor¬ 
responding highest maximum mass for pure quark stars 
in Fig. [12] with g^j = 2 and = 300 MeV held fixed 
is located at Wo- ^ 725 MeV and ^ 26 MeV with 
a value of ^ I.SMq which is not compatible with the 
recent pulsar mass measurements. 

A high mass of the sigma meson seems to be necessary 
in order to fulfill the constraints for pure quark star con¬ 
figurations, i.e. equation (l40l) and (jlll), in contrast to the 
variation of the vector coupling constant. A projection 
on the (7(^-TOo--plane on the other hand (with a fixed value 
of B^/"^ = 120 MeV) leads to a null result for pure quark 
star configurations due to the high value of the vacuum 
constant. Only a small value of the vacuum constant 
leads to pure quark star configurations. 

As the two solar mass limit for pure quark stars with 
g^j = 2.0 is not reached, we increase the repulsive vector 
interaction g^j in the following to a value of g^j = 3.0. 
The resulting maximum masses and lines of constraints 
are shown in Fig. [T31 There is an overall increase of the 
maximum mass as expected for a greater repulsive inter¬ 
action. Surprisingly, for rria > 600 MeV the maximum 
mass increases with the a meson mass again, which can 
not be seen in Fig. [T^l This behavior indicates a switch 
in the dominance of the scalar- and vector field contribu¬ 
tions to the stiffness of the equation of state. Above a cer¬ 
tain value, around g^; ^ 2.75, the repulsive fields gain on 
their influence on the maximum masses compared to the 
attractive scalar fields. A combination of g^j ^ 2.75 and 
rria > 600 MeV leads then to higher maximum masses 
with increasing a meson mass instead. 

The lines of the two constraints intersect at B^^^ ^ 
32 MeV and nia ^ 760 MeV. From this point on pure 
quark star configurations are realized between the two 
lines of constraint as discussed above. Note that in 
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Fig. [13] results are shown for cr meson masses of up to 
= 1000 MeV. At the intersecting point at ~ 
32 MeV and ra^j ~ 760 MeV the maximum mass of 
the quark star would be ^ 2.7Mq, being smaller for 
larger values of the tr meson mass Wo- and larger val¬ 
ues of the vacuum constant. A pure quark star with a 
maximum mass of ~ 2.0Mq can be found at a vacuum 
constant of = 70 MeV for a meson masses between 
900 MeV <m„ < 1000 MeV. 

F. Comparison with other models 

In the following we compare our findings with other 
approaches for studying pure quark stars, i.e. selfbound 
strange stars. For hybrid stars not discussed here a low- 
density hadronic equation of state needs to be augmented 
for a thorough discussion which is beyond the scope of 
the present investigation, but will be adressed in |52l| . 

Coelho et al. use a SU(2)-flavor symmetric NJL 
model with a repulsive vector coupling. For a large vector 
coupling their EoS stiffens like in our model calculations 
leading to higher masses at given radii. 

Weissenborn et al. use an extended quark bag 
model. Strange stars can reach maximum masses be¬ 
yond 2Mq in their work if additional terms compared 
to the standard MIT bag model are introduced, either 
in the form of some effective interaction motivated from 
one-gluon exchange or from a large gap motivated from 
color-superconductivity. They found a maximum mass 
for a pure quark star to be at 2.7>AMq, which is in the 
same order of magnitude as in this work. 

In the work by Torres and Menezes [s^ pure quark 
stars would have maximum masses of 2.29Mq. They use 
a model, where the quark masses are assumed to have a 
certain given density dependence causing a stiffening of 
the EoS compared to the standard MIT bag model. 

V. SUMMARY 

In this work a chiral Quark-Meson model in SU(3)- 
flavor symmetry has been studied for the description 
of compact stars consisting of pure quark matter. The 
thermodynamical properties have been calculated via the 
grand potential in the zero temperature limit. The gap 
equations were solved selfconsistently to determine the 
EoS, which serves as an input to solve the TOV equations 


and compute the mass-radius relations. The EoS and the 
dependence on the four free parameters of the model, the 
vector coupling constant, the constituent quark mass, the 
sigma meson, and the vacuum constant, were systemati¬ 
cally investigated. 

The variation of the vector coupling constant showed 
the highest impact on the EoS. The higher its value, the 
stiffer the EoS, leading to maximum masses in excess of 
2Mq. The dependence of the EoS with the constituent 
quark mass mg in vacuum, which fixes the scalar coupling 
constant, is such that the smaller mg, the stiffer the EoS. 
The mass of the a-meson studied covers a range from 
400 to 1000 MeV. Eor a smaller mass of the tr-meson the 
EoS becomes stiffer. Finally, the vacuum constant does 
not affect the values of the meson fields, it just shifts the 
energy density at a given pressure. The EoS substantially 
stiffens when decreasing the vacuum constant so that for 
small values of < 40 MeV maximum masses of up 
to > 2Mq could easily be reached. 

The stability of two-flavor and three-flavor quark mat¬ 
ter have been checked as well, to see whether or not the 
model parameter space allows for physically reasonable 
quark matter properties in the SU(3)-flavor approach. 
The 2Mq limit set by the recently discovered millisec¬ 
ond pulsars PSR J1614-2230 [3 and PSR J0348-b0432 
[3 could be reached. Having considered the stability 
constraints in the equations (HOI) and (HD, most choices 
of the parameter space were found to be hybrid stars. 
Nonetheless, pure quark star configurations with ~ 2Mq 
can be realized in a small physically reasonable param¬ 
eter range. A sizable additional repulsion mediated by 
the exchange of vector mesons as well as a nonvanishing 
vacuum pressure seems to be crucial to allow for maxi¬ 
mum mass configurations of quark stars compatible with 
the recent pulsar mass measurements. 
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